Unifying quantum heat transfer and superradiant signature in a nonequilibrium collective-qubit system: A polaron-transformed Redfield approach
1. IntroductionUnderstanding the mechanism of nonequilibrium quantum transport in low dimensional systems is a long-standing challenge, which has been extensively investigated from solid-state physics,[1] quantum thermodynamics,[2] molecular electronics[3] to quantum biology.[4] In particular, quantum thermal transport, where the particle and heat flows are modulated by the temperature bias, has triggered the emergence of phonics.[5–7] Phonons have been successfully utilized to fabricate various functional devices, such as thermal diode, memory, and transistor.[8–11] In analogy to phonics, quantum heat transfer mediated by the spin (e.g., anharmonic molecule or qubit) has also been intensively analyzed in a similar way, which leads to the realization of quantum spin-thermal transistor,[12,13] spin heat engine,[14–16] and nonequilibrium spin network.[17,18]
As a generic model to describe quantum heat transfer at the nanoscale, the nonequilibrium spin-boson model (NESB) is composed of a two-level-system (TLS) that is coupled to two thermal baths,[19] which was originally proposed to study the quantum dissipation.[20–22] Many approaches have been explored to investigate the underlying mechanism of heat transfer in the NESB.[23–34] Analytically, the Redfield scheme is properly introduced to investigate the sequential process in the weak spin–bath interaction regime, where two thermal baths show additive contributions to the heat transfer.[14,35] However, the Redfield approach breaks down in the strong spin–bath interaction regime, where the multiphonon excitations should be involved to characterize the nonequilibrium heat-exchange. Then, the nonequilibrium noninteracting blip approximation (NIBA) can be applied to study baths induced nonadditive and cooperative transfer processes.[23,36–39] However, both methods are found to have their own limitations; i.e., the Redfield approach is unable to capture the multi-phonon processes in the strong coupling regime, and the heat current based on the nonequilibrium NIBA scheme does not show linear proportion to the coupling strength in the weak coupling regime.[38] Recently, the nonequilibrium polaron-transformed Redfield equation (NE-PTRE) was proposed to successfully unify the steady state heat current in the NESB.[32,40,41] However, an exploration of the NE-PTRE to study more nonequilibrium spin systems is lacking but is urgently required for the spin-based quantum heat transfer. This paper aims to fill this gap by applying the NE-PTRE to analyze the heat transfer in the collective-qubit model.
Recently, a superradiant signature of quantum heat transfer has been discovered in the collective-qubit system with weak qubit–bath interaction.[42,43] The steady state superradiance describes the effect that qubits collectively exchange energy with thermal baths, resulting in the current scaling as
, with Ns the number of qubits. In sharp contrast, the superradiant heat transfer vanishes in the strong coupling regime based on the nonequilibrium NIBA under the Marcus approximation.[44] Hence, the steady state behavior of the collective-qubit system is significant distinct from each other in limiting interaction regimes. It is consequently demanding to analyze the heat transfer feature from weak to strong couplings from a unified perspective.
To address these problems, we extend the NE-PTRE combined with full counting statistics (FCS)[45,46] to investigate quantum heat transfer in the collective-qubit system. The counting-field dependent NE-PTRE successfully unifies the current and fluctuations (e.g., heat current, current noise, skewness), with the incoherent picture in the weak quit–bath coupling limit and the multi-phonon excited transfer picture in the strong coupling limit. Moreover, superradiant heat transfer is investigated at large temperature bias,[42] and the disappearance of the superradiant signature is explained by enlarging the number of qubits beyond the weak qubit–bath interaction.
This paper is organized as follows. In Section 2, the collective-qubit model and the NE-PTRE are described. In Section 3, FCS is briefly introduced and the counting-field dependent NE-PTRE is derived, which enables us to analyze the steady state heat transfer. In Section 4, the heat current, current noise, and skewness are all found to hold the unified features by extending the application of the NE-PTRE. In Section 5, we study the transition of superradiant heat transfer from weak to strong couplings, and explain the vanishing mechanism of superradiant signature with large number of qubits. In the final section, we present a concise summary.
2. Nonequilibrium collective-qubit system2.1. ModelThe nonequilibrium energy transfer in the collective-qubit system, which interacts with two thermal baths, is modeled as
. The collective-qubit model is described as
where the collective angular-momentum operators are
with
being the Pauli operator of the
i-th qubit,
Ns is the number of qubits,
ɛ and Δ are the Zeeman splitting energy and the coherent tunneling strength of the angular momentum operator, respectively. In the limit of
Ns=1, the Hamiltonian in Eq. (
1) is reduced to the seminal nonequilibrium spin-boson model.
[19,35] The
u-th thermal bath is composed of noninteracting bosons, shown as
, where
creates (annihilates) one phonon with frequency
. The system–bath interaction is given by
where
is the coupling strength between the angular-momentum and the
u-th thermal bath, and is characterized by the spectral function
. In this paper, we select the spectral function having the super-Ohmic form
, where
is the coupling strength and
is the cutoff frequency of the thermal baths. The super-Ohmic spectrum has been extensively considered to investigate the quantum dissipative dynamics and transport in molecular electronics,
[47–49] solid-state devices,
[50] and light-harvesting complexes.
[51,52]To analyze the multi-phonon involved energy transfer processes, we apply the canonical transformation
to the Hamiltonian
as
,[53,54] where the collective phonon momentum is
. After the transformation, the modified system Hamiltonian is given by
where the factor
is specified as
and the renormalization energy is
. In weak system–bath interaction limit
, the renormalization factor becomes
η=1 and
. While in strong system–bath interaction region
, the factor
η approaches zero, and the renormalization energy
. The generation of large
ξ will dramatically change the energy structure of the modified qubits system, compared to the counterpart in the weak coupling limit. Moreover, the modified system–bath interaction is given by
2.2. Nonequilibrium polaron-transformed Redfield equationWe apply the NE-PTRE to investigate the dynamics of the collective-qubit model. The NE-PTRE, which is one type of the quantum master equation, has been successfully applied to unify the nonequilibrium energy transfer in the seminal spin-boson model.[32,41] It is known that the modified system–bath interaction disappears under the thermal average (
). Hence, it may be safe to perturb
up to the second order to obtain the quantum master equation. Under the Born approximation, the density operator of the whole system can be decomposed as
, where
is the density operator of the qubits and
is the density operator of the baths at equilibrium. The quantum master equation in the Markovian limit is obtained as
where the correlation functions are
with the phonon propagator in
u-th thermal bath
. Furthermore, in the eigen-basis
, the dynamical equation can be re-expressed as
where the transition rate is
and the element is
. The rate
describes that even (odd) number phonons are involved in the transfer process between the states
and
.
In the weak qubit–bath coupling limit, the heat transfer is dominated by the sequential process and
. Thus, the correlation function
is reduced to
, and
by ignoring the high-order correlations. Accordingly, the lowest order of the transition rate
includes the term
, with
. Moreover, all of the off-diagonal elements of the qubits density matrix approach zero at the steady state (not shown here, which is quite similar to the result in Fig. 1). Considering the commutating relation
, the quantum master equation in Eq. (8) after long time evolution is simplified as
which is identical to the dynamical equation based on the Redfield scheme (see Appendix
A for the details).
In the strong qubit–bath interaction limit, the coherent tunneling of qubits in Eq. (3) is dramatically suppressed (
), and the correlation factor
vanishes. However, the other factor
is kept finite, which contributes to the quantum heat transfer. Hence, the correlation functions in Eq. (7) are reduced to
. Consequently, the master equation in Eq. (8) on the local basis
with
is changed to
where the rate is
, with
, the energy gap
, and the correlation function in the frequency domain
. Hence, the dynamical equation is completely reduced to the nonequilibrium NIBA scheme (see Eq. (
B4) in Appendix
B for further details).
3. Full counting statistics of heat current3.1. The general theoryFull counting statistics is a two-time projection protocol to measure the current and fluctuations.[45,46] For the energy transfer in the multi-terminal setup, the generating function is generally given by[55]
where
is the counting-field parameter to count the energy flow into the
u-th bath with the Hamiltonian
,
with the propagator
, and
is the initial density matrix of the whole system. Moreover, considering the modified propagator
, it can be re-expressed as
, with
. Hence, the generating function is re-expressed as
After the long-time evolution, the cumulant generating function is obtained as
Therefore, the
n-th cumulant of heat current fluctuations into the
u-th bath is given by
Specifically, the lowest cumulant is the steady state heat current
, the second cumulant is the zero-frequency current noise
, and the third cumulant is the current skewness
. In the following, we apply the theory of the FCS to obtain the counting-field dependent NE-PTRE of the nonequilibrium collective-qubit model.
3.2. Counting-field dependent NE-PTRETo count the heat flow into the right-hand thermal bath, we introduce the counting field parameters as
,[55] which results in
where the system–bath interaction with the counting field parameter is expressed as
with
and
. Then, after the canonical transformation
with transformation operator
and
the modified Hamiltonian is given by
where the modified system–bath interaction with counting field parameter is given by
By perturbing
under the Born–Markov approximation, we obtain the second-order quantum master equation as
where the correlation functions with counting field parameter are
In absence of the counting field parameter, the correlation functions are reduced to the standard version in Eq. (
7). In the eigen-basis, the dynamics of the density matrix elements can be specified as
where the transition rates are
If we reorganize the reduced density matrix of collective-qubit from the matrix form to the vector expression, the dynamical equation in Eq. (21) is expressed as
where
is the superoperator to dominate the system evolution. Then, the cumulant generating function at
t-time is given by
where
is the vector form of the initial system density matrix, and
is the left-eigenvector of
as
, with the normalization relation
. Hence, heat current fluctuations at the steady state can be straightforwardly obtained by following Eq. (
14).
4. Unified steady state heat currentQuantum heat transfer in the NESB has been successfully investigated based on the Redfield and noninteracting-blip approximation schemes in the weak and strong qubit–bath coupling limits, respectively. However, the steady state heat current was found to be distinct from each other in a broad coupling regime.[38] Recently, the nonequilibrium polaron-transformed Redfield equation combined with full counting statistics was proposed to unify the heat current between these two limiting approaches.[32,41]
Here, we try to extend the counting-field dependent NE-PTRE to unify the heat transfer in the nonequilibrium collective-qubit model in Fig. 1. At resonance (ɛ=0), we first analyze the steady state heat current in Figs. 1(a) and 1(d). It clearly exhibits the turnover behavior, which unifies the counterparts within the Redfield and NIBA methods as the qubit–bath coupling strength approaches the weak and strong coupling limits. Although we admit that to gain an analytical expression of the heat current at arbitrary qubit–bath coupling is rather difficult, it can be obtained in limiting regimes. Here, we study the analytical expression of the steady state heat current with Ns=2. Specifically, in the weak interaction limit, the heat current is analytically expressed as (see Eq. (A9) in Appendix A)
where the coefficient
which is linearly proportional to the qubit–bath coupling strength. The current in Eq. (
25) is the special case of the general expression in Eq. (
A9). It is found in Fig.
1(a) that
becomes identical to the counterpart from NE-PTRE in the weak coupling limit (e.g.,
α=0.001). Moreover, it should be noted that the current
Jweak in Eq. (
25) with
Ns=2 has a similar expression to the case in the standard NESB (
Ns=1) in the weak coupling limit, which are both proportional to the thermodynamic bias (i.e.,
).
[23]While in the strong coupling limit, the dynamical equation in Eq. (19) with the number of qubits Ns=2 is reduced to the kinetic form (see Eq. (B4) in Appendix B)
where the population is
with
. The transition rate is
where the coefficient is
, the energy gap
, and the correlation function in the frequency domain is
with the factor
If we select
Ns=1 (
), then the rate
is reduced to the standard NESB result (see Eq. (
20) in Ref. [
36]). Thus, the cumulant generating function at steady state is given by
where
. Consequently, the heat current is given by
where the first (second) term describes the process that as the qubits release (gain) energy
ξ, the right-hand bath absorbs (emits) phonon energy
ω and the left-hand bath obtains (provides) the remaining energy
.
Jstrong shows a similar structure with the counterpart in the standard NESB, which is jointly contributed by the two thermal baths.
[36] Moreover, the expression of
Jstrong captures the turnover feature of heat current in Figs.
1(a) and
1(d) as shown by the dashed-blue lines with squares. Hence, we conclude that the steady state heat flux is clearly unified in the nonequilibrium collective-qubit model.
Next, we analyze the zero-frequency current noise and skewness in Figs. 1(b), 1(c), 1(e), and 1(f). It is interesting to find that the results based on the NE-PTRE also perfectly bridge the limiting counterparts in the weak and strong coupling regimes, which may demonstrate the unification of current fluctuations in the extended spin-boson systems (e.g., collective-qubit model). Moreover, the turnover behavior which is presented in the current is unraveled for the second and the third cumulants. Although not shown here, higher cumulants of current fluctuations also show such unified features. We should note that all above results are valid both at resonant and off-resonant conditions, which clearly exhibits that full counting statistics of the heat current is generally unified within the NE-PTRE scheme.
5. Suppression of superradiant heat transferThe superradiant effect, which describes that the system exhibits collective response under the modulation of the external field or thermal bath, has been extensively investigated in quantum phase transition,[56,57] critical heat engine,[58] and quantum transport.[42,43] In particular, quantum heat transfer in the nonequilibrium large-spin system shows the superradiant signature in the weak spin–bath coupling regime.[42] Specifically, under a large temperature bias (e.g.,
), the steady state heat current is expressed as
which follows the condition
. It should be noted that equation (
30) is a special case of the current in Eq. (
A9) in Appendix
A. However, with the strong spin–bath interaction based on the NIBA scheme combined with the Marcus approximation, it is found that such superradiant transfer vanishes.
[44] Hence, we apply the NE-PTRE to clarify this apparent paradox.
We first study the effect of the temperature bias
on the heat current in Fig. 2(a) with spin-bath coupling strength
, which is considered weak for the seminal spin-boson model (
). It is found that the current shows monotonic enhancement by increasing both
and Ns. Moreover, at large temperature bias (e.g., TL=8 and TR=0.4), the current becomes nearly stable with large number of qubits, which is clearly exhibited in Fig. 2(b) (e.g., Ns=32). While for relatively small number of qubits (e.g.,
) in Fig. 2(b), the superradiant signature of heat current is numerically obtained as
with
(for both resonance and off-resonance). Thus, it is known that the expression of superradiant heat current in Eq. (30) becomes invalid at large Ns.
To exploit the origin, we plot the current as a function of the qubit–bath coupling strength α for different Ns with TL=8, TR=0.4, and ɛ=0, as shown in Fig. 2(c). It is shown that for small Ns (e.g., Ns=8), the heat flux exhibits approximately linear increasing behavior in the weak coupling regime (e.g.,
). The heat flux demonstrates the sequential transfer process, where a superradiant heat transfer is accordingly observed, which could be described by the Redfield scheme. While the current shows distinct behavior for large Ns from that for small Ns in the same coupling regime. For example, the current for Ns=32 with the coupling strength α=0.01 has already surpassed the turnover point of the current, while the current for Ns=8 is almost linearly dependent on the coupling strength in the same coupling regime. It is known that the appearance of the turnover point of the current is the significant signature of the multi-phonons involved coherent transfer, as exploited in the nonequilibrium spin-boson model.[32] Since η is nearly equal to 1 in these two cases and the heat current shows different features in the same coupling regime, we conclude that
does not necessarily correspond to the weak coupling condition. Furthermore, the behavior of the current should be properly described by the NIBA scheme, which results in the absence of nonequilibrium superradiant signature. Therefore, the effect of superradiant heat transfer in the collective-qubit model will be dramatically suppressed in the large Ns regime; i.e. it is an Ns-dependent phenomenon.
6. ConclusionIn summary, we investigate the quantum heat transfer in the nonequilibrium collective-qubit system by applying the nonequilibrium polaron-transformed Redfield equation combined with full counting statistics. We first analyze the effect of qubit–bath coupling on the steady state heat current, which results in a turnover behavior and is similar to the counterpart in the nonequilibrium spin-boson model. Interestingly, the current consistently bridges the results in the weak and strong coupling limits, which clearly demonstrates that the heat current can be unified in the multi-qubits case. Although it should be admitted that the general solution of the heat current is difficult to obtain even for Ns=2, the analytical expression can be still obtained in the weak and strong couplings based on the Redfield (Eq. (25)) and NIBA (Eq. (29)) schemes, respectively. Moreover, the current noise and skewness are shown to be unified accordingly. Although not shown in the present paper, the unification of higher cumulants of heat current can also be observed. Therefore, we propose that full counting statistics of heat current at steady state can be unified by using the NE-PTRE.
Next, we study the superradiant heat transfer in the high temperature bias regime. It is found that with small number of qubits, the heat transfer is described by the sequential process under the Redfield scheme. The heat current exhibits
, an apparent signature of the steady state superradiance. While with the large number of qubits, the superradiant signature of the heat flux vanishes. The corresponding physical process is described by the NIBA scheme because multi-phonons should be involved to contribute to the heat transfer. Therefore, we conclude that the superradiant transfer in the collective-qubit model is a size-dependent phenomenon, and it will be strongly suppressed by increasing the qubits number.
We believe that based on the counting-field dependent NE-PTRE, the unified feature of steady state heat transfer may be realized in a much bigger family of the quantum nonequilibrium system, such as a quantum spin-boson network.[18]